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Quantum states are inherently fragile, making their storage a major concern for many practical 
applications and experimental tests of quantum mechanics. The field of quantum memories is 
concerned with how this storage may be achieved, covering everything from the physical systems 
best suited to the task to the abstract methods that may be used to increase performance. This 
review concerns itself with the latter, giving an overview of error correction and self-correction, and 
how they may be used to achieve fault-tolerant quantum computation. The planar code is presented 
as a concrete example, both as a quantum memory and as a framework for quantum computation. 



INTRODUCTION 

What is a quantum memory? The answer to this ques- 
tion is straightforward: it is something that can be used 
to store a quantum state. But a quantum state of what? 
Stored for what purpose? These questions lead to a 
multiplicity of answers. And this multiplicity leads to 
many diverse fields of research throughout quantum op- 
tics, quantum information and condensed matter [JHH]. 

Research into quantum memories can be split into two 
main fields: the physical and the logical. In reality the 
split is not so neat with some approaches belonging to 
both, but it is a useful fiction for the purpose of this 
review. The field of physical quantum memories concerns 
the actual physical systems that may be used to store a 
quantum state, systems such as quantum dots, rare-earth 
ions in solids, NV centres in diamond, Raman memories, 
cold trapped atomic ensembles and many others. In- 
depth reviews of such systems can be found in P^HTT] . 

Each of these approaches has its own strengths and 
weaknesses. They allow different operations to be per- 
formed on the stored information, and different ways to 
interface the memory with other systems. The physical 
system used for a quantum memory will therefore depend 
on the application that the memory is to be used for. 

Some applications, such as certain proposals for quan- 
tum computation, require quantum memories to store 
states with very high fidelity for long times. These high 
requirements can be beyond what we can hope any phys- 
ical memory to achieve on its own. It will need some 
help to keep errors in check, to keep the fidelity high and 
elongate the decoherence time. It's research into how this 
help may be achieved that is the focus of research into log- 
ical quantum memories. The typical approach is that of 
multipartite encoding. Instead of storing a single qubit 
in a single physical memory, a large number of physi- 
cal memories are used. The qubit, now called a logical 
qubit, is then stored non-locally across the entire multi- 
memory system. Any error on one memory can be seen 
by comparing its state to that of its neighbours. This is 
done either through measurement, allowing the error to 
be removed by quantum error correction |19| . or by im- 
plementing a Hamiltonian to energetically suppress the 



errors. Many proposals have been put forward for such 
logical encodings, capable of removing multiple forms of 
errors including spurious operations on systems as well 
as loss of some physical memories. A (by no means com- 
plete) list of some examples of error correction in simple 
systems can be found here [20H24] . 

The restrictions and requirements of the physical mem- 
ories will affect what sort of encoding is possible for the 
logical qubits, as well as what sort of encoding is required 
to best protect against the typical errors of the system. 
As such, it is important for there to be a deep under- 
standing between those who work in these fields. It is a 
study of these the logical quantum memories, and their 
place in quantum computation, that concerns the bulk 
of this review. In particular we consider the capabilities 
of Kiteav's surface codes, which include the famous toric 
code, to act as both a quantum memory, and as a fault- 
tolerant platform for quantum computation |1H I12j . 

WHY SHOULD WE STORE A QUANTUM 
STATE? 

First let us consider a few of the applications of a quan- 
tum memory, with particular emphasis on those applica- 
tions that require the use of a logical encoding on top of 
the physical one. 

Quantum memories and flying qubits 

Many applications of quantum information theory re- 
quire the use of so-called flying qubits. These typi- 
cally take the form of photons, with the qubit state 
encoded within the photon polarization. Flying qubits 
allow quantum information to be transferred over po- 
tentially long distances, providing the basic resource for 
many important applications such as quantum communi- 
cation protocols, experimental probes of quantum foun- 
dations and the transfer of quantum information between 
different parts of a quantum computer. 

Even so, it isn't enough just to have photons whizzing 
around. The qubits need to be processed in order to 
extract useful information, and undo the effects of errors 
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that occur during their journeys. For this, they need to 
be caught, stored, manipulated, measured and released 
in a deterministic manner. These are the tasks for which 
quantum memories are required. The exact relationship 
that quantum memories need to have with flying qubits 
depends on the application. We will now consider a few 
examples. 

Quantum repeaters 

Many proposed applications of quantum technology 
rely on the possibility of long-range quantum state trans- 
fer, including quantum key distribution [TH3| and quan- 
tum networks [HI US] ■ Just as long-range transmission of 
classical information is crucial in present day communica- 
tions, the equivalent quantum process is likely to become 
the workhorse of the future. However, transmission of 
quantum information over long distances will always in- 
cur heavy losses. And, unlike in the classical case, copy- 
ing of the information to amplify the signal is forbidden 
by the no-cloning theorem [53 1^ . As such, a quantum 
solution must be provided for this quantum problem. 

A solution, proposed by Hans Briegel et al. in 1998, is 
known as the quantum repeater approach For this, 
a number of quantum repeater stations are set up be- 
tween the sender and recipient. The long distance over 
which the information is to be sent is therefore split up 
into many smaller links. Each station will then create 
entangled pairs, keeping one particle of each and sending 
the other to one of its neighbours. By doing so, neigh- 
bouring stations will share entanglement. Entanglement 
swapping [37| [55] and purification [5S] may then be used 
to extend the range of this entanglement, ultimately pro- 
viding a number of entangled pairs shared by the sender 
and recipient. Of course, the noise incurred when par- 
ticles are sent between neighbouring stations will mean 
that this entanglement is imperfect. Nevertheless, the 
use of purification can distil many noisy copies into one 
of high fidelity. With this, a quantum state may be tele- 
ported directly from sender to recipient without error 
[50] . even though no particle has had to travel further 
than the distance between two stations. 

Quantum memories are crucial to the implementation 
of a quantum repeater. Without storage of the entan- 
gled states created and received by each station, the en- 
tanglement between them cannot be maintained. It also 
cannot be manipulated to create the purified entangle- 
ment that stretches the entire distance of the network, 
required for the ultimate transmission by teleportation. 
For these things to be achieved, it must be possible to 
store the quantum states of the photons sent between 
neighbouring stations. The quantum memory used for 
this must also have a long coherence time, and allow the 
operations required for entanglement swapping and state 
purification. To achieve this, quantum repeaters may 



employ quantum error correction [31H33j . 

Quantum foundations 

The weight of experimental evidence supports the fact 
that Bell's inequalities are violated in nature [Ml [55] . 
Nevertheless, all current experiments have suffered from 
loopholes that still allow interpretations consistent with 
local realism to be made. It is therefore of fundamental 
interest to close the gaps as much as possible, to further 
demonstrate the non-local nature of entanglement |36| . 
Quantum memories help toward this end in two ways. 
Firstly they aid the preparation of high fidelity entan- 
gled pairs separated by the long distances required to 
ensure locality. This can be done, for example, using the 
quantum repeater technique described above. Secondly 
they provide more reliable interfacing with the quantum 
state, allowing measurements to be made with the high 
efficiencies required to close the fair sampling loophole. 

Another cornerstone of quantum mechanics is the un- 
certainty principle. However, an entropic reformulation 
of this has led to a more detailed understanding of exactly 
what can and cannot be concurrently known. Specifi- 
cally, it is found that entangling a system to a quantum 
memory allows measurements to be made that are ap- 
parently more accurate than Heisenberg allows [37]. Of 
course this is not the case, and the results will always 
be in perfect agreement with the generalized uncertainty 
relationship that takes the entanglement into account, 
but it demonstrates the power of quantum memories in 
probing the dark corners of quantum mechanics. 

Quantum computation 

Quantum computation aims to quickly perform com- 
putations for which no efficient algorithm classical com- 
putation algorithm is known, and for which it may be 
that none exist. Through the preparation, manipulation 
and processing of quantum information, it should be pos- 
sible to efficiently solve such problems as quantum sim- 
ulation [38] and the factoring of primes [39], as well as 
allow significant speed-ups for other tasks [^ . Just as in 
classical computers, the ability to store quantum infor- 
mation during the computation, to manipulate it while 
stored and to send it elsewhere when required, will be a 
critical element of any architecture. Also, such memories 
must be of very high fidelity. Otherwise, the build up 
of errors during the computation will yield a useless re- 
sult. The development of a quantum RAM [3T1I32], and 
quantum hard drives for long term storage, is therefore 
a major focus of research into quantum computing. 

Because of the many varied and stringent conditions 
for the storage of information within quantum comput- 
ers, the use of error correction to create logical quantum 
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memories is required [15] • It may, in fact, be advanta- 
geous for the logical qubit states to be initialized, pro- 
cessed and read out, all without ever leaving the error 
correcting code [i5Hi5] . 

QUANTUM ERROR CORRECTION 

So far we have spoken of the possibility of combining 
many physical memories into a logical memory of higher 
fidelity. But how may this be achieved? In any quantum 
memory we typically have access to a system with a large 
number of degrees of freedom, and we wish to store a rel- 
atively small (typically two-level) quantum state. The 
approach is then to find a two dimensional subspace of 
the system for which the effects of decoherence is par- 
ticularly small, and then use this to encode the state we 
wish to store. By a clever choice of both the system and 
the subspace, we therefore allow our state to survive a lot 
longer than if it were encoded in a more exposed manner. 

It is exactly this approach that is used with error cor- 
recting codes, in both the classical and quantum realm 
[T^. To demonstrate the idea behind them, we start 
with the simplest example of classical error correction: 
the repetition code. 

Let us suppose that we have a logical bit that we wish 
to store, and the physical means to store single bit values, 

and 1. However, this physical memory is noisy. By the 
time we come to reading out the information, there is a 
probability p that a bit flip has occurred. The most naive 
way to store our logical bit is simply to encode it within 
a single physical memory. However, this method allows 
the noise to directly affect the logical bit, flipping it with 
probability p and destroying the integrity of the stored 
information. 

The easiest way to suppress the effects of the noise 
is simply to copy the logical bit three times. So, if the 
logical bit in state 0, we prepare three physical memories 
in state 000, and for logical bit 1 we prepare 111. At the 
time of readout, we may find that the three physical bits 
are in the same state, 000 for example. The probability 
that they started in state 0, and no flips occurred, is 
{l—pY- The probability that they started in state 1, and 
a flip occurred on each, is p^ . Clearly (assuming p < 0.5), 
the former is more likely than the latter. We may also 
find two of the physical bits in the same state and one in 
the other, 001 for example. The probability that they all 
started in state and suffered a single flip is clearly more 
likely than the probability that they all started in state 

1 and suffered two flips. This allows us to do a majority 
voting to correct errors. If the majority of the physical 
bits are in a certain state, we assume they all started in 
that state and hence that was the state of the logical bit. 
The probability that this process fails is related to the 
smallest error that can turn one majority into the other, 
i.e. two bit flips. The probability of a logical error is 



therefore reduced from p, to ^ p^ . If we used more than 
three physical bits per logical bit, the error rate would be 
reduced further. In fact, the logical state can be retrieved 
with arbitrary fidelity for any p < 0.5 as the number of 
physical bits per logical bit is increased. 

The repetition code is the classical worlds natural 
means of storing classical information. It is simply an 
abstract expression of the fact that classical information 
is stored in states of many particles, such as the many 
molecules of ink used to print a letter in a book. Even if 
the letter is smudged or parts of it are obscured, it has 
to be perturbed significantly before a passage cannot be 
read. 

The repetition code has a natural extension to qubit 
encoding only if the quantum state is known. It is only 
in this case that multiple qubits can be prepared in the 
same state. However, this is a trivial case, since a known 
quantum state can be stored easily in a classical computer 
simply by recording each entry of the density matrix. In 
any practical application of quantum information, it is 
the storage of arbitrary quantum states without any prior 
knowledge that is important. In this case the no-cloning 
theorem applies [251 126j , preventing such straightforward 
repetition encoding. 

One possible quantum generalization could then be as 
follows. Say we have an arbitrary qubit state | ip) — 
a I 0) + /3 1 1) (a mixed state would also be valid, but we 
use pure for simplicity). We can then take two ancilla 
qubits, each in state | 0), and apply a CNOT with each 
ancilla as target and the original qubit as source. The 
resulting state is a \ 000) + /3 | 111) , and so | 0) is encoded 
in I 000) and 1 1) in |111). The repetition encoding in 
the computational basis means that a physical bit flip 
error rate of p becomes a logical bit flip rate of ~ p^. 
However, a phase flip acting on any qubit will cause a 
phase flip to the logical qubit, and no measurable trace 
of the error will be left. A physical phase flip error rate of 
p therefore becomes a logical phase flip rate of ^ 3p. The 
encoding protects against bit flips, but allows phase flips 
to occur more easily. Obviously this does not provide a 
good quantum code. 

A solution to this problem is to concatenate the repe- 
tition code in two orthogonal bases. To explain this, we 
will use the plus/minus basis, the eigenstates of the Pauli 
X operator on qubits, 

I±> = ^(|0)±|1)). (1) 

We will first use three qubits to encode logical plus/minus 
states I +)3 = I + + +) and | — )3 = | ) with a stan- 
dard repetition encoding. The computational basis states 
in this encoding will take the form, 

|0)3 = ^(l + + +) + l )) (2) 

1 1)3 = 4^(1 + + +)-! ))• (3) 



4 



This encoding is just that described above, but with the 
roles of bit and phase flips interchanged. Here it is the 
phase flips that are suppressed, with a logical rate of 
^ and the bit flips whose rate is increased to ^ 3p. 
To suppress the latter we use three groups of three qubits 
to give a repetition encoding of the computational basis, 

10)9 = 10)310)310)3, 11)9 = 11)311)311)3- (4) 
The total encoding is then, 

1 0) 9 = ^(1 + + +) + ! ))^^ (5) 

1 1) 9 = ^(1 +++)-! — (6) 

Both bit and phase flips are suppressed with error rates 
going from p ^ p^, leading to a lower logical error rate 
than the physical error rate. Since bit and phase flips are 
the Pauli operators X and Z, respectively (making Y a 
conjunction of the two) and since all single qubit errors 
can be decomposed into single qubit Pauli's, this code 
therefore suppresses all errors. Such suppression in all 
bases makes this a good quantum error correcting code. 
It is one of the simplest to achieve such suppression, and 
is known as the 9-qubit code [46], for obvious reasons. 

From this example, we see that classical error correct- 
ing codes do not work for quantum information. How- 
ever, mashing classical codes together in two orthogonal 
bases can work. This provides a powerful method of cre- 
ating quantum codes [HI |47l [48] , but are there quantum 
codes that do not take this form? Can we find a more 
general framework? 

To begin a more general definition of quantum error 
correcting codes, let us recall the motivation mentioned 
above: we wish to take a large system, find a subspace 
for which the effects of decoherence is particularly small, 
and then use this to encode a logical qubit. To do this 
we will take a large number of physical quantum memo- 
ries, whose joint Hilbert space will have the large number 
of degrees of freedom required. We then choose a sub- 
space of states spread non-locally throughout the phys- 
ical memories. This is chosen such that, if the state of 
the system is moved out of the subspace by any small 
error, the signature of this error can be detected by a 
cleverly defined set of measurements without disturbing 
the stored state. These are known as syndrome mea- 
surements. Once these have revealed that the error has 
occurred, it can either be undone or its effects can be 
accounted for in any subsequent measurements. 

An important class of quantum error correcting codes 
are the stabilizer codes [TOl 149) . To define these, one first 
takes a number of physical memories which each encode 
2-level systems, and then chooses a number of mutually 
commuting operators acting on these. These operators 
are known as the stabilizers of the code, they are all 
products of Pauli operators and so all have eigenvalues 



-1-1 and —1. The stabilizer space of the code is defined to 
be the mutual -1-1 eigenspace of all stabilizers. It is this 
subspace that is used to encode the logical state. 

Note that the stabilizer operators are Hermitian, and 
hence are observables. We may therefore consider a mea- 
surement of all stabilizers, with each giving a result of 
-f 1 and —1. Since the state of the system should always 
reside in the stabilizer space, the results of all these mea- 
surements should be -1-1. The presence of a stabilizer with 
outcome — 1 is therefore a signature that an error has oc- 
curred, and so the stabilizer measurements play the role 
of syndrome measurements. For a well defined stabilizer 
code, any local error will disturb some stabilizers in this 
way. An analysis of the results of a measurement of all 
stabilizers, known as the syndrome, should then be able 
to undo the effects of the errors by returning the state not 
only to the stabilizer space, but to its initial state within 
this space before the error occurred. The resilience of the 
logical state to decoherence is then much better than the 
states of the physical memories. 

In the next section we will consider a concrete example 
of a stabilizer code: the planar code. However, it should 
be noted that there are many approaches to quantum er- 
ror correction other than the stabilizer formalism. Some 
are equivalent to stabilizers or generalizations of the ap- 
proach, and some not. These include CSS codes [47ll48]. 
qudit stabilizer codes [5DH5^ , stabilizer codes with weak- 
ened conditions [53], subsystem codes and topological 
subsystem codes [5iH57] , decoherence free subspaces [SS] 
and autonomous error correction [59H62] . 

PLANAR CODE QUANTUM COMPUTATION 

One particularly interesting form of stabilizer codes are 
topological stabilizer codes, such as Kitaev's famous toric 
code [m [T^. For these the physical memories are ar- 
ranged on a two-dimensional surface, with stabilizers de- 
fined on local clusters of these surrounding various points 
throughout the surface. The eigenstates of the stabilizers 
are interpreted as occupancies of so-called anyonic quasi- 
particles, with the eigenvalue -1-1 denoting no anyon is 
present on the region for which the stabilizer is defined, 
and —1 denoting that an anyon is present. The stabilizer 
space is therefore the anyonic vacuum. The effect of local 
errors is to locally create pairs of anyons, and the stabi- 
lizer measurement is a measurement of whether anyons 
exist at each point. Correction is achieved by trying to 
determine in which pairs the anyons were created (up to 
an equivalence class) and reannihilating them. 

We will now consider the planar code as one particular 
example of a topological stabilizer code. This will be ex- 
plained first at an abstract level in terms of the anyonic 
quasiparticles. At this level, we will see how quantum 
information can be stored and even how it may be ma- 
nipulated for universal quantum computation. We will 
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FIG. 1. The exchange statistics of (a) bosons and (b) 
fermions. 



then move to more concrete ground, describing how it 
may be reahzed on a spin lattice. This allows us to look 
at the effects of errors, the means for their suppression 
and correction and review of the error thresholds that 
have been found. Some readers may prefer to review 
this concrete section before starting on the abstract. Fi- 
nally we will review the idea of using the planar code as 
a self-correcting quantum memory, by implementing an 
additional Hamiltonian. 

It should be noted that there are multiple different 
ways in which the planar code may be used for quantum 
computation. Including transversal proposals |12| and 
ones that involve cutting and gluing of codes jSHj- The 
primary family of proposals, however, are those that use 
and manipulate 'hole'-like defects in the code. These 
were first introduced by Raussendorf, et al., [43H45] . 
though similar ideas existed previously [51]. It is the 
Raussedorf type approach that we will focus on in this 
review. However, for reasons of clarity, we will focus on 
a single-hole encoding of qubits, rather than the original 
two-hole encoding. 



Bosons, fermions and anyons 

The properties of bosons and fermions arc familiar to 
any physicist. Their exchange behaviour is summarized 
in Fig. [T] If one has a state of two bosons, and then ex- 
changes them with an exchange operator R, the resultant 
state will be exactly the same as the initial one. This is 
due to the indistinguishable nature of bosons and the fact 
that the exchange results a trivial phase (i.e., it results 
in a phase of 1). For fermions the situation is the same, 
except that the exchange leads to a phase of —1. 

Two successive exchanges of two particles is equivalent 
to a full monodromy, where one of the particles is braided 
fully around the other. In both the bosonic and fermionic 
case the single exchange phase squares to unity. The 
effects of a full monodromy is therefore trivial in both 
cases. This is a requirement placed on any point particles 



FIG. 2. The exchange statistics of e and m anyons, with e 
anyons shown in red and m in blue. 



in a universe with three spacial dimensions, since a loop 
around a point in 3D space can be continuously deformed 
to a loop that is not around a point. A full monodromy is 
therefore topologically equivalent to a trivial operation, 
and must have a trivial phase. 

In a universe of two spacial dimensions, this restriction 
no longer holds. A loop around a point can no longer be 
continuously deformed to a loop that is not around a 
point, and hence a full monodromy of point particles can 
have a non-trivial effect. Such particles would be neither 
bosons or fermions, and so we call them anyons. 



An abstract overview of the planar code 

In the planar code we will imagine that we have access 
to a 2D 'universe' in which there are only two elemen- 
tary particles, which we call e and m. The e is its own 
antiparticle, and accumulates a trivial phase when ex- 
changed with another e. The same is true for the m's. 
However, their anyonic nature becomes evident when we 
consider their exchanges with each other. A full mon- 
odromy of an e around an m (or equivalently, vice-versa) 
yields a phase of —1, as shown in Fig. [2j 

We will consider this universe to be square in shape, 
as depicted in Fig. [3] The anyons e and m, as point 
particles, can exist at any localized point in the universe. 
However, it is also possible for them to exist in a delocal- 
ized way on the edges of the universe. Specifically, the 
top edge of the universe is able to hold a single e anyon. 
This anyon, if present, cannot be thought of as being at 
any specific point along the edge. Even if one sees it move 
from a localized point in the bulk and 'fall' of the edge, 
one cannot think of it as being just beyond the edge from 
where it disappeared. Rather, the anyonic occupation is 
spread out over the entire edge. An e anyon can simi- 
larly exist on the bottom edge, and an m anyon on each 
of the two sides. The top and bottom edges appear as 
hard boundaries to m anyons, as do the sides to e's. An 
example of a configuration of anyons in this universe can 
be seen in Fig. [3j 

As we will see in the next section, in reality this uni- 
verse will be a spin lattice plus a little imagination. We 
will define operators acting on the spins of the lattice, 
which will correspond to the stabilizers of an error cor- 
recting code. The properties of these operators will allow 



FIG. 3. The 'universe' on which the e and m anyons hve, 
with an example configuration of anyons on the bulk. In 
general, a configuration of any number of anyons on the bulk is 
possible. The occupations of the edges will always ensure that 
the total number of both e and m anyons is even. Note that 
the lattice background to this figure and the two which follow 
is not important in itself. It is only to show the distinction 
between the 'universe' and the empty space that surrounds 
it, and to foreshadow our later discussion in terms of a spin 
lattice. 



us to interpret their measurement outcomes as occupa- 
tions of quasiparticles, with an outcome of —1 meaning 
a particle is present in the region on which the operator 
is defined, and +1 meaning that one is not. The two 
particle types stem from two different types of stabilizer. 
By performing operations on the spins of the lattice, we 
gain certain powers over the anyons of our universe. At 
this level of abstraction, these powers may seem quite ar- 
bitrary. However, when we consider what is really hap- 
pening in terms of the underlying spin lattice, we will 
see how they may be justified. The first power we allow 
ourselves is the ability to create pairs of anyons from the 
vacuum. We can create a pair of e or to anyons (but 
not an e and an to since they are not antiparticles) , and 
also any arbitrary superposition of creating such pairs 
and creating nothing. We also allow ourselves to move 
anyons, throw them off their respective edges and anni- 
hilate pairs of anyons. We can also measure the anyon 
occupancy of any point on the bulk as well as the edges, 
and hence know where all the anyons in the universe are 
whenever we wish. 

The final powers we allow ourselves concern 'holes' in 
the universe, as in Fig. |4] There are two types of holes 
that can exist, which we call rough and smooth, due to 
their appearance when the spin lattice realization of the 



FIG. 4. The universe with four holes, two rough (left and 
right) and two smooth (top and bottom) and an example 
configuration of anyons. 



model is considered. We allow ourselves to create both 
types of hole and alter their shape as we wish, allowing 
us also to move them around. These holes behave in a 
similar way as edges. Rough (smooth) holes can accom- 
modate a single delocalized e (to) anyon, but appear as 
hard boundaries to m (e) anyons. As such, we also al- 
low ourselves to measure the anyonic occupations of the 
holes. 

So, now we have such a universe, is there anything in- 
teresting that we can do with it? It turns out that there 
is: we can use it to perform quantum computation. The 
first step towards this is to first determine a means of 
storing the quantum information that we wish to com- 
pute. As usual we will do this in the form of qubits: 
two level systems with states labelled | 0) and 1 1) in the 
so-called computational basis and | -f ) and | — ) in the or- 
thogonal plus/minus basis. The most basic unitary op- 
erations that can be performed on these are the logical 
Z and X operations, Pauli operators which act on the 
basis states in the following manner, 

^|0) = |0), Z|l) = -|l), 
X|+) = |+), X\-) = -\-), 
X|0) = |1>, X|1> = |0), 
Z|+) = |->, Z\-) = \+). 

We will use holes to store qubits. For each smooth hole 
we store a single qubit (a smooth qubit), associating the 
two possible anyon occupations, vacuum and to anyon, 
with the logical qubit states | 0) and 1 1), respectively. For 
each rough hole we also store a qubit (a rough qubit). 
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associating vacuum and e anyon occupations with the 
states I +) and | — ), respectively. In order to satisfy the 
restriction to even numbers of both e and m anyons, any 
spare anyons are moved off the edge. The bulk is kept 
anyon free. 

The logical X operation maps a qubit in state | 0) to 
one in state | 1), and vice- versa. To implement this on a 
smooth qubit, we therefore need to map a hole contain- 
ing no anyon to one containing an anyon, and vice- versa. 
The corresponding operation will obviously be one which 
creates a pair of m anyons, moves one into the hole, and 
moves the other off the edge. If the hole previously con- 
tained vacuum, it will contain an anyon after this process. 
If it contained an anyon before, the process annihilates 
this and results in the vacuum. Similarly, a logical Z 
can be performed on a rough qubit by creating an e pair, 
placing one in the hole and moving the other off the edge. 

Such operations need to act along all parts of the uni- 
verse in a path that stretches from the hole to an edge. 
It is this that gives the qubit its topological protection 
against errors. No small, local error can effect the qubit 
as long as the hole is kept far from the edge. Such an 
error could create a pair of m anyons on the edge of 
the hole, and cause one to fall in. This would change 
the anyonic occupation of the hole, and hence appear to 
have caused an X error. However, the presence of the 
other anyon near the hole is a clear signature of what 
happened. Any measurement of the anyon configuration 
would reveal that the error occurred, and the spare anyon 
may be thrown in the hole to annihilate its partner and 
undo the problem. 

Logical Z operations for rough qubits and logical X 
operations for smooth qubits require a phase of — 1 to be 
generated if an anyon resides in the associated hole, and 
no phase to be generated otherwise. To achieve this, we 
can use the braiding properties of the anyons. To apply 
a, Z to a smooth qubit we first create a pair of e anyons, 
braid one around the hole, and then annihilate the pair. 
If the hole contained vacuum (and hence the qubit was 
in state | 0), this process has trivial effect. If it contained 
an m anyon (qubit in state 1 1) ) , the braiding yields the 
required phase of —1. Similarly an X can be applied to 
rough qubits by creating an m pair, braiding one and 
annihilating. 

Such operations also need to act along a path, this time 
one which completely encircles a hole. Making the holes 
large therefore topologically protects the qubits. Any 
local error can only create and move anyons around a 
small portion of the hole. The probability of many such 
errors conspiring to braid an anyon all (or most of) the 
way around is suppressed with the hole's diameter. The 
encoding and the logical operators are depicted in Fig. 

El 

Since any logical error can be decomposed into Z and 
X errors, the fact that both of these are topologically 
protected for both qubit types (rough and smooth) is 





FIG. 5. Four holes in which four qubits (two rough and 
two small) are stored, (a) The paths along which m anyons 
must be moved in order to perform logical X operations. For 
smooth qubits these stretch from hole to edge, changing the 
occupancy of both. For rough qubits these encircle the hole, 
(b) Similarly, the paths required for e anyons to perform log- 
ical Z operations. 



sufhcient to show that the qubits are completely topo- 
logically protected so long as the holes are large and well 
separated [SSHSTj . However, a full treatment of 

errors is deferred until later. 

It should also be noted that the logical encoding allows 
fault-tolerant preparation of smooth qubits in the | 0) 
and 1 1) basis states, and rough qubits in the | and 
I — ) basis states. This is done by first clearing all anyons 
away before defining the holes. It is then known that the 
holes must contain the vacuum, and any local errors that 
may alter this will leave the same signature as described 
above. This gives the preparation of smooth | 0) states 
and rough | -|-) states, and allows preparation of smooth 
I I) states and rough | — ) states by application of an X 
and Z, respectively. 

An entangling operation can be performed on the 
stored qubits in a very straightforward manner. Since 
we are allowed to move the holes, let us consider braid- 
ing a smooth hole in a full monodromy around a rough 
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hole (or vice- versa). The effects of this obviously de- 
pend on the anyonic occupations of the two holes. If the 
smooth hole contains the vacuum (and hence is in logical 
state I 0)), the braiding will be trivial. If it contains an 
m anyon (and hence is in logical state 1 1)), the effect is 
to braid an m around the rough hole, generating a phase 
of —1 iff that hole contains an e anyon (and hence is in 
state |— )). This corresponds to a logical X operation 
on the rough qubit. The total effect of the braiding of 
the two holes is therefore to implement the entangling 
CNOT gate, with the smooth qubit as the control and 
the rough as the target. Like the logical X and Z oper- 
ators, the CNOT is also topologically protected. Local 
errors cannot cause the braiding of the hole to take place 
in error if the holes are sufficiently large. Also, though 
the operations used to move the holes may be imperfect 
and cause local creation and movement of anyons, these 
can be identified and corrected as long as the movement 
is sufficiently slow. 

The ability to realize a fault-tolerant entangling gate 
demonstrates the power that this model has to process 
quantum information. However, the gate set that we can 
realize through the creation and braiding of the anyons is 
still not universal. For that we require something more. 

To see what we can use to achieve universality, con- 
sider in general some non- universal gate set, such as the 
Clifford group of operators. Add to this set another op- 
erator, namely preparation of a given single qubit state. 
If this new set of operators is universal, the state can be 
called a 'magic state' [Til [72] ■ The abihty to prepare this 
state and use it as an ancilla boosts a non-universal gate 
set to universal. 

Given the gates that we can already achieve, it can be 
shown that the required magic states for smooth qubits 
are of the form, 

\9) = ^{\O)+expz0\l)), (7) 

for = 0, 7r/4, 7r/2, tt. With similar states required for 
rough qubits, except that | 0) is changed for | and | 1) 
for I — ). To prepare these, smooth qubits can be moved 
near to side edges, shortening the distance required to 
perform the logical X operator. In fact, using the spin 
lattice realization of the next section, it will allow the 
logical X basis to be accessed through operations on a 
single spin |65l l68H70j . The manipulations required to 
prepare a state such as the above may then be achieved 
through X-basis measurements. Some of these may lead 
to the creation of unwanted e particles, but these will be 
removed through error correction. 

Of course, making something easy for us to do means 
making it easy for nature also. If we can perform oper- 
ations in the X basis, so too can nature perform errors. 
However, note that the states with = 0, tt are eigen- 
states of the X operator, and hence are immune to X er- 
rors. We needn't worry about these as long as the holes 



are large enough to prevent Z errors. The same does 
not apply to the 6 = 7r/4, 7r/2 states. The errors can 
be minimized by moving the holes away from the edge 
as soon as the preparation is finished, but a significant 
amount will still occur. To deal with these, it is found 
that distillation can be used ptiiSl [7T1 [7^ . This takes 
many noisy copies of the states and uses the original (and 
non-universal) gates available to the model to distil these 
down to a single copy with arbitrary accuracy. 

The magic states can therefore be prepared fault- 
tolerantly. Added to the other fault-tolerant gates possi- 
ble through the manipulation of the anyons, this allows 
the model to achieve the holy grail of fault-tolerant uni- 
versal quantum computation. 

Realizing the planar code on a spin lattice 

In order for the ideas developed above to bear fruit, 
we must find some way to simulate this anyonic universe 
in the lab. To do this, we need a lattice of two-level 
quantum systems. It is usual to call these systems 'spins', 
and hence we use a 'spin lattice'. However, the systems 
need not actually have any spin associated with them. 
Any two-level quantum systems will do, so long as they 
can be manipulated in the required way. As such, it is in 
fact photon based systems that are responsible for most 
experimental progress in the planar code so far |73fl75j . 

It is here that the physical quantum memories enter the 
discussion, since it is these that will serve as the spins. To 
implement the planar code we take many physical quan- 
tum memories, each of which have a noise rate greater 
than the required level, and place them in a lattice such 
as that of Fig. [6] The size of this lattice is characterized 
by the number of plaquettes along the horizontal direc- 
tion or the number of vertices along the vertical. These 
should be of the same order. We will take them to be 
equal, and denote them by L. 

Multi-spin observables are defined on the spins, and 
correspond to occupations of e and m on different points 
on the lattice. The steps developed in the previous sec- 
tion can then be used to create a number of logical quan- 
tum memories, in which qubits can be stored at a greatly 
reduced error rate and used for fault-tolerant quantum 
computation. 

The two bases we will use when studying the states 
of the spins are the computational basis and plus/minus 
basis, as used for the logical qubits above. The Pauli 
operators for the spins will be denoted cr^, and cr^. 

The observables used to define the occupancies of m 
anyons are associated with each plaquette (or face) of 
the lattice. Specifically, consider the spins of the lattice 
in the computational basis. Each plaquette is surrounded 
by four spins (three for those near the edge). If, for a 
given plaquette, an even number of these spins are in the 
state 1 1), we say that that plaquette does not hold an 
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FIG. 6. The lattice on which the planar code is defined. A 
selection of single spin errors and their efTects in creating and 
moving anyons are also shown. Note that the edge anyons 
denote anyons that have been moved off the edge by errors, 
rather than the edge occupations themselves. 



anyon. If there are an odd number of | l)'s, we say that 
it holds an m anyon. 

Similarly, the occupancies of e anyons are defined on 
the vertices of the lattice, which are also surrounded by 
four spins, or three near the edge. If, for a given vertex, 
an even number of these four spins are in the state | — ), 
we say that that vertex does not hold an anyon. If there 
are an odd number of | — )'s, we say that it holds an e 
anyon. 

The occupations of the edges are defined in the same 
way. If an even number of the spins on the left (right) 
edge are the state the left (right) edge holds no 
anyons. If an odd number, it holds an m anyon. An 
even number of | — )'s on the top (bottom) edge means no 
anyon and an odd number means an e anyon is present. 

The definition of the anyons is no more complex or 
mysterious than this. The e and m anyons simply cor- 
respond to the parity of | — )'s and | l)'s around vertices 
and plaquettes, respectively. From these simple defini- 
tions, we find that the anyons will have all the properties 
described in the previous section. However, before we 
continue, let us state the definitions in a more formal 
way. To do that, the following observables are defined 
for the spins around each plaquette, p, and vertex, s, 

As^\{af, B, = \{at, (8) 

where i € s denotes the spins i that surround the vertex 
s, etc. These are the stabilizer operators of the code. 
Each observable Bp has eigenvalues -1-1 and —1. The +1 
(— 1) eigenspace is formed by states for which there is an 



even (odd) number of | l)'s around the plaquette p. The 
-1-1 and —1 eigenspaces therefore correspond to states of 
vacuum and an m anyon on the plaquette p, respectively. 
These eigenvalues for the Ag similarly correspond to the 
vacuum and an e anyon on the vertex s, respectively. 

Note that these operators all mutually commute. This 
allows the configuration of e and m anyons to both be well 
defined at the same time. Measurement of e occupations 
on vertices will not disturb m occupations on plaquettes, 
or vice-versa. The stabilizer space of the code is defined 
as the space of states that satisfy the conditions, 

As\ij) = \^P}, i?p I V') = I V-) . (9) 

The stabilizer space is therefore the mutual +1 
eigenspace of all the plaquette and vertex operators, and 
corresponds to the anyonic vacuum. 

In addition to the stabilizer operators are the edge op- 
erators, which act on the stabilizer space. The left and 
right edge operators are products of ct^'s down the left 
and right edges of the code, respectively. The top and 
bottom operators are products of 's along the top and 
bottom. These commute with all stabilizers, and the 
product of left and right commutes with the product of 
top and bottom to ensure that the net anyon occupation 
off the entire edge is well defined, however they do not 
mutually commute. This will be important for the edge 
encoding considered later. 

Let us consider the system in a state within the stabi- 
lizer space, and look at the effects of local errors. Since 
all operators can be decomposed into the Pauli basis, it 
is sufficient to consider only Pauli operators. A ap- 
plied to a spin i will anticommute with two Bp operators, 
those for the plaquettes that neighbour the spin. The ef- 
fect is therefore to flip the sign of these two stabilizers, 
creating an m anyon on the two plaquettes. The oper- 
ator commutes with all other stabilizers, and so has no 
further effect. The application of a cr^ on any spin will 
similarly create e anyons on neighbouring vertices, and a 
will create both an e and an m pair. 

Now let us take a state that is not within the stabi- 
lizer space, but already has an m anyon present on the 
plaquette p. What is the effect of a applied to the 
spin shared by the plaquette p and its neighbour p'7 It 
will flip the sign of both, returning the state to the -1-1 
eigenspace of Bp, but placing it in the —1 eigenspace of 
Bpi. This cr^ therefore moves the m anyon, and u^'s will 
similarly move e anyons. 

In general, the effects of Pauli operators are to create, 
move and annihilate anyons, as shown in Fig. [6] Note 
that to create a pair of m anyons, move them and then 
annihilate them, a loop of cTx operations is required. This 
loop will be equal to a product of the As stabilizers for 
the vertices enclosed by a loop. The same is true for e 
anyons, except that the loop will be one of Bp stabilizers. 
Since the stabilizer space is an eigenstate of all stabilizers, 
such closed anyon loops will have trivial effect. 
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The anyonic braiding of the e and m anyons is due 
entirely to the anticommutation of cr^ and ax operations. 
For example, consider the creation and movement of an 
e pair by the application of a string of a^s. Creating a 
m pair, braiding one of them around of the e's, and then 
reannihilating the m pair, can be achieved by a loop of 
operators that is equivalent to a product of Ag stabilizers. 
The fact of the braiding means that the loop will enclose 
the e that was braided around, and so the loop of cr^'s 
will overlap with the string of cr^'s on an odd number 
of spins. The resulting anticommutation results in the 
phase of — 1 acquired through the braiding. 

To store qubits in the way described in our earlier ab- 
straction we need to make holes in the code. If the holes 
were to be static, this would be done simply by modifying 
the lattice on which the code is defined. Instead of the 
standard square lattice, a square lattices with holes could 
be used. However, since we wish holes to be a resource 
that can be created and moved at will, we will need to 
specify a means by which this can be done. 

Let us begin with the standard square lattice planar 
code. Creating a hole means choosing an area of the 
lattice on which we no longer wish to enforce stabilizers. 
This will mean that the spins within that region will not 
be part of the code for as long as the hole resides in that 
place. As such, they must be disentangled from the rest 
of the system. 

To create a smooth hole, we first choose a selection of 
plaquette stabilizers that we no longer wish to enforce, 
and then determine which spins are acted upon only by 
those plaquettes and no others. These are then disentan- 
gled by performing measurements in the basis, and 
hence the plus/minus basis, of the spins. To keep things 
tidy, all spins are rotated (if needed) to the | state 
after this measurement. 

There will be, in general, some vertex stabilizers with 
support only on the disentangled spins. These will be 
of vertices that reside within the hole. These stabilizers 
are also no longer enforced. There will also be stabilizers 
with partial support on the disentangled spins. Specifi- 
cally, those residing on all points of the perimeter of the 
hole, but not on the corners, will act on one disentangled 
spin. The definition of these stabilizers is then changed 
accordingly, from four-body operators to three body op- 
erators acting only on the spins in the code. If the original 
stabilizer was in its -|-1 eigenspace, and hence the vertex 
held no anyon, the same may not be true for the reduced 
stabilizer. The creation of the smooth hole will therefore 
result in spare e anyons that need to be cleared up. The 
end result is then to create a large plaquette in the code, 
which is the composite of the smaller plaquettes whose 
stabilizers are no longer enforced. Like any other pla- 
quette, this one can either hold vacuum or an m anyon. 
However, unlike the others, it has no stabilizer associated 
with it. As such, the dimension of the stabilizer space is 
increased by the presence of the hole, since it holds states 




FIG. 7. The creation and movement of a smooth hole is 
shown, first by creation, then expansion and finally contrac- 
tion. Creation is performed by measuring the spins within 
the hole in the basis. The parity of the number of spins 
with outcome | — ) will be equal to the parity of the number 
of e's initially within the vertices contained by the hole, and 
hence will be even if the state was initially in the stabilizer 
space. For every spin with outcome | — ), the corresponding 
truncated stabilizer will hold an e anyon. These must be re- 
moved, or otherwise taken into account, and the spins in state 
I — ) should be rotated to | -|-). Expansion is simply achieved 
through more creation. The additional e anyons will only arise 
on the expanded boundary. To contract a hole, previously un- 
enforced stabilizers need to be measured. If these contain m 
anyons, they must be placed in the hole. If they contain e 
anyons, these must be annihilated. Creation and movement 
of rough holes are achieved similarly, with the roles of plaque- 
ttes and vertices, and of e and m anyons, reversed. Also, the 
basis is used in place of the . 



associated with both possible occupations. The creation 
of a smooth hole is depicted in Fig. [7] 

If the code was initially in the stabilizer space, with 
anyonic vacuum everywhere, the hole will be initialized in 
the vacuum state. There will always be an even number 
of spare e anyons created, which can be annihilated with 
each other without affecting the hole state. In general, 
the hole will hold vacuum if the component plaquettes 
initially held net vacuum (an even number of m's) and 
will hold an m if the component plaquettes initially held 
a net m. 

The creation of a rough hole is similar. A selection 
of vertices are chosen for which the stabilizers will no 
longer be enforced. The spins touched only by these are 
identified and measured in the computational basis of 
the spins. Stabilizers for plaquettes within the holes are 
no longer enforced, and those bordering it are reduced if 
required. This will result in an even number of spare m 
anyons to be cleared up. The hole is then effectively a 
large vertex (though it is not often depicted as such for 
sake of clarity) that is the composite of those vertices that 
are no longer enforced. It can hold either the vacuum 
or an e anyon, with its initial state dependent on the 
previous occupations of the component vertices. 

Since the presence of a hole increases the stabilizer 
space, they can be used to store quantum information. 
Specifically, each hole can store a logical qubit. The | 0) 
and 1 1) states of a smooth hole are defined by the vacuum 
and m occupancies, respectively. The logical Z operator 
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FIG. 8. Logical operators for smooth holes. The Z operator 
is a loop of 's around the hole, and the X is a string of 's 
from the hole to the edge. For rough holes the roles of Z and 
X, and and , are interchanged. 



is therefore the product of the (now no longer enforced) 
stabilizers of the component plaquettes, and takes the 
form of a product of az operators that encircle the spin. 
Similarly the logical X for rough holes, where the | +) 
and I — ) states are defined by the presence or lack of an 
e anyon, is the product of stabilizers for the component 
vertices. 

Note that it was earlier stated that products of stabi- 
lizers act trivially on the stabilizer space, and yet here 
we find that such a product is a logical operator. The 
resolution to this apparent contradiction is that the sta- 
bilizers for the component plaquettes and vertices are no 
longer a part of the code. It is only products of stabi- 
lizers within the code that act trivially. Anyon loops are 
therefore not trivial if the anyons have braided around 
holes, but are trivial otherwise. 

The logical X operator for smooth holes requires the 
occupancy of the hole to be changed. It does this by 
creating a pair of m anyons, and placing one within the 
hole. However, the requirement for logical states to lie 
within the stabilizer space means that the other m anyon 
cannot reside on any plaquette for which the stabilizer is 
still enforced. The logical operator must therefore take 
this spare anyon to, and place it over, the edge of the 
code. The logical X is therefore a string of operators 
from a smooth hole to the left or right edge. The logical Z 
for rough holes is similarly a string of ct^'s from the rough 
hole to the top or bottom edge. These logical operators 
are depicted in Fig. [Sj 

The movement of a hole is achieved through its expan- 
sion and contraction. To expand a smooth hole, simply 
choose more plaquette stabilizers on the boundary of the 
existing hole that you wish to no longer enforce. Iden- 
tify the spins that are no part of no enforced plaquette 
stabilizer, and measure them in the basis. Finally re- 
define the reduced vertex stabilizers on the edge of the 
hole and measure them, finding the spare e anyons that 
have arisen on these and removing them. The movement 
of a smooth hole is depicted in Fig. [7] 

For this final step, there is an additional complication 



than must be considered. When a hole is being moved, 
its logical qubit will typically be in an arbitrary state. As 
such, the application of unwanted logical Z errors must 
be avoided. However, consider the case that the hole is 
expanded in all directions, and so retains none of its ini- 
tial boundary. Spare e anyons will be created randomly 
around the entire perimeter of the hole. Two completely 
valid ways to annihilate these with each other will be 
possible. However, the operators required to implement 
these differ by a factor of the Z operator. As such, we 
are equally likely to annihilate the anyons the right way, 
which does not implement a Z on the stored qubit, and 
the wrong way which does. 

On the other hand, consider an expansion that only 
pushes back part of the hole's perimeter. In this case it 
is known that the operator for annihilation of the spare 
anyons should only have support on that part of the 
perimeter that was expanded. The choice of how to anni- 
hilate is then unique (up to trivial factors of stabilizers) 
and will not implement a logical Z. As such, it must 
be remembered that holes are to be expanded slowly. A 
small portion must be expanded, then error correction 
must be run to annihilate any anyons that have arisen 
during the expansion, either because of the expansion or 
because of spin errors, then the hole can be expanded a 
little more. All the same arguments also apply to rough 
holes. 



Edge encoding 

The hole based encoding is required to perform the 
quantum computation scheme described in the earlier ab- 
straction. However, there is a simpler encoding that can 
be used if the planar code is required only to serve as a 
quantum memory. It's simplicity also means that it is 
the basis of proof-of-principle experiments. 

The simplest way to store qubits in this code is using 
the edge occupations. A logical | 0) can be stored in the 
state with vacuum on the left and right edges, and | 1) in 
the state with an m anyon on these edges. The logical 
I +) (I — )) is stored using vacuum (e anyons) on the top 
and bottom edges. This is depicted in Fig. [9j Because 
the occupation operator for the top or bottom edges does 
not commute with that for the left or right, these are not 
independent qubits as they would be for a hole based en- 
coding. Instead they are the logical Z and X basis states 
of a single qubit. The logical X and Z operators corre- 
spond to the edge occupation operators and so stretch 
across the entire code. In order to perform a logical er- 
ror, nature must therefore create a pair of anyons and 
transport them the width of the code, with each going 
to opposite edges. As for the hole based encoding, these 
means that errors are suppressed for large system size. 

Though the hole based encoding is the primary ap- 
proach for quantum computation in the planar code. 
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FIG. 9. A depiction of the computational basis encoding 
using the edges. The state | 0) is defined as a state with no 
m anyons on the left and right edges, and | 1) as a state with 
m anyons on these edges. The top and bottom edges are 
correspondingly a superposition of both holding vacuum and 
both holding e. 



other approaches have been proposed using this edge en- 
coding. These include transversal application of gates 
[T^ . and so-called 'lattice surgery', which interacts the 
logical qubits stored on different codes by combining the 
codes and then splitting them |63j . 



Stabilizer measurements 

By now it should be obvious that the workhorse of 
the planar code, when used solely as a quantum memory 
or when quantum computation is performed, is measure- 
ment. It is by measurement that errors are detected such 
that they may be corrected, and it is measurement that 
makes and moves the holes. All other operations need 
only act on single spins. It is the measurements that 
drive their effects non-locally and allow them to have 
computational power. 

Most important of the measurements are the stabi- 
lizer measurements, which are the only multi-spin aspect 
of the proposal. Since there are four-spin observables, 
their measurement is obviously not straightforward. The 
naive approach would be to deduce the outcome of the 
four-spin measurement by making four single-spin mea- 
surements. Measurement of for example, determines 



whether there are an odd or even number of | l)'s around 
the plaquette p. This could be deduced by measuring 
each of these spins in the computational basis and count- 
ing the number of | l)'s. However, this approach would 
not commute with the neighbouring operators, and 
so would cause errors. Measuring all plaquette sta- 
bilizers in this way would then leads to logical Z errors, 
as would measuring all vertex stabilizers like this lead to 
logical X errors. 

A means to measure the stabilizers must be therefore 
be used that is both realistic and commutes with all other 
stabilizers. For the most part, direct measurement of the 
four-spin operator is not realistic, though some excep- 
tions exist jl06| I118j . The generally favoured approach 
is therefore to use ancilla spins located at the centre of 
every vertex and plaquette. Code spins then couple to 
these using CNOT gates in such a way that the outcome 
of the four-spin measurement is mapped onto the single 
spin ancilla state. Single spin measurement may then be 
used to reveal its value. 

To be more specific, an anyon is located on a plaquette 
if the state is in the —1 eigenspace of Bp, and so if there 
are an odd number of spins around the plaquette in state 
I 1). To determine this, we take a spin initialized in state 
I 0) and place it inside the plaquette. CNOT's may then 
be performed with each spin around the plaquette as the 
source and the ancilla as the target. A CNOT performs a 
CT^ on the target only if the source is in state | 1) , and so 
as many 's will be performed on the ancilla as there are 
I l)'s around the plaquette. An even number will cancel 
each other out, leaving the ancilla in state | 0). An odd 
number acts as a single cr^, leaving the ancilla in state 
I 1). Measuring the ancilla and obtaining the result | 0) 
or 1 1) then means that the plaquette holds the vacuum 
or an anyon, respectively. 

For vertices, it is the parity of | — )'s that determines 
whether an anyon is present. As such, a phase controlled 
NOT is required, which performs a on the target only 
if the source is in state | — ). This is obtained from a 
CNOT simply by performing a Hadamard on the source 
spin before and after the CNOT's application. The ver- 
tex occupation is read-out similarly to the above by plac- 
ing a I 0) initialized ancilla in the center of the vertex, 
and performing phase controlled NOT's with each ver- 
tex spin as the source and the ancilla of the target. The 
final ancilla state then represents whether there was vac- 
uum in the vertex, or an anyon. Note that the CNOT's 
and phase controlled NOT's do not individually commute 
with the surrounding stabilizers. However, once the four 
operations required for any plaquette or vertex are ap- 
plied, the composite operation does commute with all 
stabilizers. This form of measurement therefore does not 
disturb other plaquettes and vertices, or cause logical er- 
rors (at least, when done perfectly). 

With this in mind, we can reduce the entire planar 
code down to a few simple ingredients. We need to place 
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a large number of physical memories in a two-dimensional 
lattice. We need to be able to act on these with single 
spin operations (measurements and rotations in the cr^ 
and bases, and the Hadamard). And we need to be 
able to perform CNOT's between neighbouring physical 
memories. If these can be performed experimentally, the 
code can be implemented. This will allow us to achieve 
either a straightforward quantum memory, or quantum 
computation. This is subject of course to the noise level, 
as we will now discuss. 



Errors and thresholds 

When considering quantum memories, one important 
figure of merit to be considered is the lifetime. This will 
depend on various parameters, such as the temperature 
of the thermal bath in which the memory is immersed, 
and the strength of the coupling to the bath. Since error 
correcting codes typically do not employ a Hamiltonian, 
we can consider their to be no correlations between the 
errors caused to each spin by the bath. This allows us to 
use a more fundamental quantity in place of the lifetime: 
the threshold probability that a errors occurs on the spins 
of the code. 

When defining error thresholds for codes, we usually 
begin by assuming that the code is prepared perfectly in 
its stabilizer space. For the planar code, this means the 
anyonic vacuum. It is then coupled to a thermal bath, 
causing single spin errors to be applied randomly and 
independently to each spin with a certain rate. As an 
example, consider a bath that causes bit flips (cr^ errors) 
to occur on each spin with a rate 7. The probability that 
a bit flip occurs within an interval dt is then 7(it, and 
the probability that a net bit flip (that is an odd number 
of bit flips, since even numbers cancel) occurs within an 
interval t is p(i) = (1 - e-T*)/2. 

For small the probability that an error occurs on any 
spin will be small. The errors therefore result in the cre- 
ation of isolated pairs of anyons. If we can assume that 
the stabilizers of the code, and hence the anyon occupan- 
cies, can be measured perfectly, we can easily see where 
these pairs are and determine which spins suffered an 
error. The correction of the errors is therefore straight- 
forward. After long however, the probability of error 
will be close to 1/2. Close to half of plaquettes and ver- 
tices will therefore hold an anyon, and it will not be clear 
which spin errors occurred to cause the configuration of 
anyons. Correction then becomes impossible. For exam- 
ples, see Fig. 
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It is known that, between these two extremes, we will 
find a threshold probability Pc- If the rate of errors is 
below this, then error correction is always possible and 
the stored state can be retrieved with unit probability 
as L — > (X). For error rates above this, correction can- 
not be performed, and fidelity of the stored state will 



FIG. 10. Example anyon configurations are shown for low 
(left) and high (right) rates for bit flip errors. In both cases, 
there was initially no m anyon on either the left or right edges, 
meaning a logical | 0) was stored. For the former configura- 
tion, the pairing of anyons is simple, and it is easy to see on 
which spins the errors occurred. These are noted on the fig- 
ure. It can been seen that the errors changed the occupancy 
of the left edge, and this can be corrected , returning the state 
to logical I 0) . For the latter, the odd number of anyons in 
the bulk makes it clear that one of the edge occupancies has 
changed. However, the configuration does not make it easy 
to determine which one. If the attempt at error correction 
chooses the wrong one, the result will be to place an m on 
both edges, leaving the state in logical | 1) and hence applying 
a logical error rather than correcting it. 



become minimal as L — > 00. This threshold probability 
then translates to a threshold time, the lifetime of the 
code, until which the stored state can be retrieved re- 
liably. If correction is performed periodically, with the 
delay between correction rounds less than the lifetime, 
the quantum memory can be sustained indefinitely. 

The retrieval of a stored state requires the use of an er- 
ror correction, or decoding, algorithm. This is a classical 
algorithm that takes the results of the stabilizer mea- 
surements and determines the best guess for what errors 
occurred. The most naive decoding algorithm is brute 
force, which considers every possible sequence of errors 
consistent with the stabilizer measurements and assesses 
their likelihood. The use of brute force will ensure that 
the success rate of the decoding is the best that is pos- 
sible, and that the correction succeeds all the way up 
to the threshold error rate. Practically, however, brute 
force is highly inefficient. The time taken to run the com- 
putation scales as expO(L^), and will take longer than 
the age of the universe even for small codes. Alternative 
algorithms are therefore required, which correct up to a 
lower threshold error rate than the theoretical maximum 
but do so with an efficient run-time. For any error cor- 
recting code, and any given error model, we will therefore 
have multiple threshold error rates. One will be the the- 
oretical maximum, obtainable by brute force, and others 
will be the thresholds obtainable by various decoding al- 
gorithms. 

The simplest error model that can be applied to the 
planar code is that of independent bit (cr^) and phase 
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(cr^) errors. Let us use p to denote the probability of 
a bit flip on each spin, and p' to denote the probabil- 
ity of a phase flip. For the planar code defined on a 
square lattice, the threshold probabilities are found to 
be Pc — p'f. = 11% [12]. Both error rates must be below 
this limit for the stored quantum state to be retrievable. 
For codes defined on other lattices, the thresholds will 
change, typically with one rising above 11% and the other 
falling below. For example, using a hexagonal lattice re- 
sults in thresholds oi p = 6.74% and p' = 16.4%. For a 
triangular lattice these two numbers will be exchanged. 
A series of randomized lattices have been found to in- 
terpolate between these extremes. In all cases it was 
found that the threshold achieves the so-called hashing 
bound, meaning that the planar code is highly optimal 
|76j . The lattice may therefore be defined to best suit 
the ratio between p and p', thus optimizing the lifetime. 
Practically, error correction may be performed using the 
minimum weight perfect matching algorithm |121 l77H79j . 
This typically achieves a threshold that falls short of the 
maximum only by around 0.5%, and so is very effective. 

Another well studied error model is that of depolariz- 
ing noise, which assumes that tr^, tr^ and tr^ errors occur 
on each spin with equal rate p/3, where p is the proba- 
bility that an error of any sort occurs. In this case, the 
threshold has been found to be Pc = 18.9% [S^. Min- 
imum weight matching is not so effective in this case, 
achieving only a threshold of 15.5%. A renormaliza- 
tion based algorithm manages 16.5%, but (theoretically) 
this can be increased by increasing the run-time [81 . A 
Markov chain Monte Carlo algorithm has been found that 
achieves an 18.5% threshold, which should be possible to 
increase to the 18.9% maximum with greater computa- 
tional cost [82] . 

The above error models assume that the stabilizers can 
be measured perfectly, giving an exact and noiseless ac- 
count of the anyon configuration created by the errors. 
However, such an idealized case is not physically relevant. 
Practically, measurement of the stabilizers will also be 
subject to errors. This means that single shot readout 
of the memory is not possible, since a single noisy mea- 
surement of the stabilizers is insufficient to perform error 
correction. Instead, the stabilizers must be measured pe- 
riodically, in order that error correction can detect not 
only the spin errors that occur, but the measurement 
errors also. 

In this spirit, an error model has been proposed and 
studied to give a realistic account of the noise that a 
system is subject to. This assumes the measurement 
scheme with CNOT's coupling to ancilla spins, as de- 
scribed above. Errors are assumed to occur on all ele- 
ments of the system. Between measurements, depolariz- 
ing noise is assumed to act on all code spins with rate 
p. The preparation of ancillary spins into their initial 
state is assumed to fail with the same probability p. The 
CNOT gates are assumed to be noisy, with single spin 



depolarizing noise acting on both spins with rate p and 
measurement is assumed to give the wrong answer with 
probability p. The current best algorithm, based on min- 
imum weight perfect matching, gives a threshold proba- 
bility oi Pa ~ 1%, the highest yet known for any local 
and two dimensional architecture [751 HH] ■ For p below 
this limit, and as L and the number of measurements are 
made large, the stored information can be read out of 
the memory with arbitrarily high fidelity and will have 
an arbitrarily long lifetime. The theoretical maximum 
for the threshold probability is not known, but devel- 
opments of the minimum weight perfect matching algo- 
rithm are constantly increasing the achieved threshold, 
and the generalization of the Markov chain Monte Carlo 
algorithm to the noisy measurement case is expected to 
yield a significant increase. 

Experiments and proposals 

The first experimental steps toward realization of the 
planar code were actually done for the toric code model, 
which is the same up to boundary conditions. In these, it 
was photons that were used for the 'spins', with entangled 
states created by parametric down conversion. The en- 
tangled states created corresponded to those of the toric 
code for small system sizes. Nevertheless, the size was 
sufhcient to allow anyons to be created and braided by 
rotations of the single photons, and for the non-trivial 
effects of the braiding to be measured [72 [71] . 

The first of these experiments realized a single- 
plaquette and four- vertex code with four photons [75] . 
A superposition of an m anyon and vacuum was created 
on the plaquette using a rotation in the tr^ basis on one 
of the spins, operations were then used to create an 
e pair on the vertices, move one of these around the pla- 
quette and reannihilate them. This changed the relative 
sign of the m-vacuum superposition in a way that could 
be measured. The second experiment followed a similar 
procedure, but with a two- vertex and four-plaquette code 
realized with six photons |74j . 

A more complex eight photon set-up has recently been 
realized using an eight-photon cluster state [75]. This 
was used to demonstrate the error-correcting properties 
of the planar code. A single qubit was encoded and shown 
to be protected from single qubit errors, with a logical 
qubit error rate lower than tthat on the physical qubits. 

Another route towards experimental realizations is 
that of cold atoms. A toolkit of methods that may be 
used to realize anyon models with optical lattices were 
explored in |83| . Experiments concerning minimal in- 
stances of topological order have also been considered 

[SHIHI]. 

One promising proposal for the full realization of the 
planar code, which would eventually allow its use as both 
a memory and a quantum computer, involves the use of 
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quantum dots [55^. The physical quantum memory in 
this case would correspond to the spin qubit of the dot. 
In order to create and manipulate the planar code states, 
some means must be found to perform the required two 
qubit gates. The traditional approach for this would be 
to use the tunnel coupling between the qubits [57] • How- 
ever, this requires the qubits to be placed very close to- 
gether, making it difficult to fit in the required gates and 
wirings. A recent proposal allows this problem to be 
avoided through the use of a long range coupling, medi- 
ated by metallic ffoating gates [55] ■ It has been shown 
that these can be used to implement CNOT gates with 
high fidelity, allowing the stabilizer measurements to be 
made on a lattice of quantum dots, and the planar code 
state to be prepared and manipulated for quantum com- 
putation. 



HAMILTONIAN PROTECTION AND 
SELF-CORRECTION 

Quantum error correcting codes are usually considered 
as a means to correct errors through measurement of 
the syndrome. However, one could also consider imple- 
menting an error-correction inspired Hamiltonian, whose 
terms refiect the syndrome measurement operators and 
whose spectrum energetically penalizes errors. Of course, 
this can only be realistically be considered for codes 
whose syndrome measurements are (quasi-)local, or such 
a Hamiltonian could never be implemented. The planar 
code (and the related toric code) is therefore a natural 
candidate for this approach, since the stabilizer operators 
of the are quasi-local, acting only on the spins clustered 
around each plaquette and vertex. A Hamiltonian built 
out of these operators may then be defined as follows, 

H = -JsY,As- JpY^Bp- (10) 

s p 

The ground state space of this Hamiltonian is the stabi- 
lizer space of the code, and the anyonic vacuum, in which 
the state of the quantum memory is stored. Anyons, 
which are created by errors, become quasiparticles whose 
creation requires an energy 2Jg for e anyons and 2Jp for 
m anyons. At zero temperature, it has been shown that 
this Hamiltonian is completely stable against (sufficiently 
weak) local perturbations. Their effects are exponentially 
suppressed with the linear size, L, of the code [89]. Th 

This is especially true when there is disorder present in 
the Hamiltonian couplings, which can induce Anderson 
localization of anyon transport [75] IUDH5^ . 

Usually the code is simply considered as a quantum 
memory when the Hamiltonian is considered, with no 
hole-based quantum computation as described above. 
However, this is not always the case. Proposals for simi- 
lar computation schemes have been made in the presence 



of the Hamiltonian, using its protection against local per- 
turbations. These are based on the planar (or toric) code 
itself [64], as well as higher dimensional generalizations 
with more complex anyonic behaviour [SSJ 1701 193] . This 
then leads into the realms of topological quantum com- 
putation with non-Abelian anyons, which is the subject 
of other reviews [MHSS]- 

At non-zero temperature the Hamiltonian greatly in- 
creases the lifetime of the code, since creation of anyon 
pairs is suppressed by a factor exp[— 0(/3 J)], for J = 
min( Js,Jp). This leads to a lifetime that scales as 
exp[0(/3J)]. The exact form of the relation for the life- 
time depends on the details of the system-bath coupling, 
the geometry of the lattice used for the code and other 
such issues. By taking such details into account, the code 
can be tailor made to ensure that the lifetime achieved 
is the highest possible |102j . 

Though this Hamiltonian gives a lifetime that increases 
exponentially with /3 and J, it is not always practical to 
raise the coupling or lower the temperature. Realistic 
regimes may therefore have a lifetime that is far below 
what we desire for quantum computation. Ideally we 
would like the lifetime to diverge with a parameter that is 
easier to increase, such as the system size. Such memories 
are called 'self-correcting' because the Hamiltonian must 
not only suppress the creation of errors, but also be able 
to undo them once they occur. 

For a classical example, we may consider the Ising 
model. For a one-dimension chain of spins of length L, 
this has the Hamiltonian, 

H.^jY^afat^,. (11) 

i 

An energy cost 2 J is given to the walls between domains 
of spin up and spin down particles, meaning the states 
with all spins up and all down are the degenerate ground 
states. Let us encode a classical bit in these states. The 
Hamiltonian is then an energetic means the enforce the 
repetition code. If coupled to a thermal bath, errors will 
cause the creation of pairs of domain walls which may 
then walk along the line, fiipping spins. 

To determine how well the Hamiltonian suppresses the 
effects of such errors, we can consider the energy barrier. 
This is the minimum amount of energy required to per- 
form a logical operation (hipping all spins in this case), 
if the operation is performed one spin at a time. The 
minimum energy required to map from the logical state 
to the logical state 1 is achieved if first a single spin is 
flipped to create two domain walls, and then neighbour- 
ing spins are flipped to move these along the line. The 
only energy cost is the initial 4J required to create the 
domain walls. The energy barrier between the two logical 
states is then 0(1) in this model. 

Consider instead the two-dimensional version of this 
model, on a square lattice with side-length L and Hamil- 
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tonian, 

H2 = jY.<yt<jt+,+a^at+L- (12) 

i 

In this case the energy cost for a domain wall is not fixed, 
but instead scales linearly with the length of the bound- 
ary. The ground states are again those with either all 
spins up or all down, and so can be used to store and 1, 
making this an alternative way to energetically enforce 
the repetition code. The main difference to the one- 
dimensional model is the mapping from to 1, which 
requires the growth of a region with opposite orientation 
until the point that it takes up the majority of the lat- 
tice. This will require it to go via a state for which the 
length of the domain wall, and hence the energy barrier, 
is 0{L). 

The different energy landscape of the two models 
causes radically different thermalization behaviour. The 
one-dimensional model has only a constant energy bar- 
rier for physical errors to become logical errors, whereas 
the two-dimensional model has a barrier which increases 
the system size. The latter is therefore far more likely 
to undo its errors than the former. This gives the one 
dimensional model a lifetime of exp[— 0(/3 J)], with no 
increase with system size. The two dimensional model, 
however, has a lifetime that increases exponentially with 
system size below the critical temperature (for the ferro- 
magnetic phase transition) and polynomially above |97j . 

Going back to the quantum regime, it is easy to see 
that the planar code Hamiltonian is a direct parallel of 
the one-dimensional Ising model. An energy is required 
to create an anyon pair, but then no further energy cost 
is required for it to perform the walk around the entire 
system that causes a logical error. This 0(1) energy cost 
of a logical error can be seen as the reason for its lifetime 
that does not increase with system size [98l |99] . For a 
self-correcting quantum memory, which we would like to 
act in a way corresponding to the two-dimensional Ising 
model, a diverging energy barrier must be constructed to 
suppress the spurious application of both logical X and 
logical Z errors. 

The most obvious way to do this is to generalize the 
planar code to higher dimensions [T^j. It has been found 
that the three-dimensional toric code has a diverging en- 
ergy cost of logical bit flip errors, but not for logical 
phase. It therefore serves as a self-correcting classical 
memory, but the lifetime of the quantum memory is 0(1) 
|100[ llOlj . The four-dimensional planar code Hamilto- 
nian does serve as a self-correcting quantum memory 
|100| . but the three-dimensional nature of our universe 
would make the experimental realization of such a model 
quite difficult. 

The greatest hope for planar code based self-correction 
therefore lies in the two-dimensional model with the ad- 
dition of long-range interactions. Multiple different forms 
of such interaction have been proposed and studied. In 



[103] a logarithmic attractive potential between anyons 
was proposed, as was a physical mechanism to realize it 
through coupling the code to an additional bosonic field. 
This leads to a self-correcting memory below a certain 
finite temperature phase transition, with lifetime scaling 
as 0{L^) with system size. Above the critical tempera- 
ture the memory will still be self-correcting, but with the 
lifetime scaling only as 0{L^) [llOj . This transition only 
occurs if the toric code is used. For the planar code, the 
lifetime will scale as 0{L^) for any temperature. 

In [76l 11041 1105j a repulsive potential between anyons 
with power law decay is considered. Though one might 
consider this repulsion to increase the rate of logical er- 
rors, it also means that the energy required to create a 
pair of anyons depends on the number of anyons already 
present. This leads to a vanishing small anyon density at 
any finite temperature. The scaling of lifetime with L is 
at most quadratic, and depends on the exact form of the 
potential. It has been found that this interaction can be 
engineered by coupling the code spins to cavity modes 
[THH] . 

In |105) a similar repulsive potential is also considered, 
but this time with an anyon- vacuum, rather than anyon- 
anyon, interaction. If the state is prepared in the anyonic 
vacuum, this means that 0(L^) vacuum plaquettes and 
vertices will oppose the creation of any anyons. This re- 
sults in an exp[— 0(L^)] suppression of anyon creation, 
and hence a lifetime that scales exp[0(i^)]. This cou- 
pling may also be realized through coupling to cavity 
modes |105) . 

Recently a novel model has been proposed, in which 
the two-dimensional toric code is coupled to a three- 
dimensional ferromagnet |lllj . All interactions in this 
model are local and of bounded strength, and no cou- 
plings to non-local features, such as cavity modes, are 
required. The effect of the ferromagnet is to induce an 
enhanced suppression of anyon creation, with the energy 
required to create a pair scaling linearly with L. This lin- 
ear energy barrier exponentially suppresses the creation 
of anyons, and so will lead to a lifetime of O(expL). 

Further proposals for self-correcting quantum memo- 
ries go beyond the planar code to completely different 
models. It has been shown that any local stabilizer 
Hamiltonian in two-dimensions jll21 I113j cannot sup- 
port self-correction, and indeed local Hamiltonians in 
any dimension will have limitations |114j . More com- 
plex approaches are therefore required. One possibility is 
that of engineered dissipation |115j . Also the Bacon-Shor 
code has been proposed as a possible candidate in three- 
dimensions [1161 1117j , however it is not known whether 
this achieves self-correction. 

Also in 3D are the cubic code |107l I108j and welded 
solid codes |109] . Both of these have the kind of energy 
barrier one would expect from a self-correcting memory. 
The minimum energy required to perform spin-by-spin 
logical operators is O(logL) in the former and 0(L^/^) 
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in the latter. However, self-correction has been shown 
not to be present for the cubic code |108j and the same is 
strongly suspected for the welded solid codes |109j . The 
problem here is one of entropy. Increasing system size 
increases the energy barrier, but also increases number 
of paths which errors may take. If the latter increases 
to fast with respect to the former, errors will be able to 
overwhelm the barrier and not self-correct. 

Even so, it can be expected that self-correcting prop- 
erties will still be present to an extent. If the system 
size is small enough with respect to the temperature, the 
number of errors present at any time will also be small, 
and so the energy barrier will have a significant effect. 
We can therefore expect the lifetime of quantum mem- 
ory stored in such codes to scale as 0{L^) up to an op- 
timal system size L* = exp[0(/3)]. For any given inverse 
temperature /?, this means a code can be constructed 
with lifetime exp[0(/3^)], by ensuring that it is of the 
right size. The lifetime does not diverge with system 
size, and so the model cannot be called self-correcting. 
However, the exp[0(/3^)] scaling of the lifetime is much 
better than the cxp[0(/3)] scaling of models with no self- 
correction properties. It is therefore called 'partially self- 
correcting'. The cubic code has been shown to behave in 
such a way |108j , opening up a new realm of possibilities 
for quantum memories. 



Experimental progress 

The first experimental realization of the toric code 
Hamiltonian was achieved using Josephson junctions 
|119| . This provided a demonstration of the degenerate 
ground state used for the quantum memory. The theory 
on which this experiment was based has a straightforward 
extension to the quantum double models [12(T . These are 
qudit generalizations of the toric code which can also be 
used to store quantum memories, and which have more 
complex anyons [llj . 

In general the major difhculty in realizing the planar 
or toric code Hamiltonian is the many-body vertex and 
plaquette interactions. Realistically, these need to be ob- 
tained perturbatively from two-body interactions. Possi- 
ble means by which this may be done are using Kitaev's 
honeycomb lattice model |121j . or perturbative gadgets 

Progress is also being made into simulations of the pla- 
nar (and toric) code model with Rydberg atoms |123j . 
Rather than realizing the Hamiltonian directly, these 
seek to simulate its dynamics using the Rydberg blockade 
|124) . Such a simulation could provide a demonstration 
of how well the quantum memory stands up to dissipative 
noise. 



Outlook 

Already we have seen a great deal of experimental 
progress in the use of error correction and other such 
protection mechanisms to enhance quantum memories. 
Many physical systems have been tamed and set to work 
storing quantum states, and small instances of quantum 
error correction have been achieved. 

Further progress, as always, requires a closing of the 
gap between theory and experiment, with concessions on 
both sides. Theory often wishes for big systems, fast 
gates and low noise, and expects measurement to be no 
harder than applying a projection operator. But none 
of these are the case, and so advances in encodings and 
protocols must be made to increase error thresholds and 
allow for the experimental complexity of performing the 
required operations. To make this easier, experiments 
must continue their push to make things bigger, faster 
and cleaner. 

For the case of self-correcting quantum memories, the 
theory is still relatively new and experiments are non- 
existent. Both matters need to be addressed. There re- 
mains theoretical work to be done in understanding what 
sorts of systems allow self-correction, and finding exam- 
ples of realistic models that could realize self-correction 
in the lab. The most advanced work so far is that done 
with long-range anyonic interactions |105L 11061 llllj . But 
it would be advantageous to have additional and inde- 
pendent approaches for comparison. A two-dimensional 
partially self-correcting memory with local interactions, 
for example, may be regarded as a favourable alternative 
despite its non-diverging lifetime with system size. In 
any case, theorists must always keep in mind that their 
definition of realistic, i.e. using only local interactions, is 
not always the same as that of experimentalists. 
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